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ABSTRACT 

Two  types  of  ionospheric  anomalies  have  been  discovered 
recently  in  the  polar  cap:  patches  and  arcs.  Polar  cap  arcs 
are  the  focus  of  this  study,  which  seeks  correlation  between 
arcs  and  total  electron  content  (TEC)  enhancements  and 
amplitude  scintillation  effects.  Simultaneous  optical  and 
radio  frequency  measurements  were  taken  at  Thule  AFB  and 
Qanaaq ,  Greenland,  using  an  All -Sky  Imaging  Photometer 
(ASIP)  and  a  specially  equipped  Global  Positioning  System 
(GPS)  receiver.  Arcs  were  discovered  to  produce 
significant,  rapidly  varying  TEC  increases,  and  small  but 
measurable  amplitude  scintillation. 
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I.   INTRODUCTION  AND  OBJECTIVES 

Current  command,  control,  communications  and 
intelligence  (C3I)  systems  operate  over  the  entire  electro- 
magnetic spectrum,  and  the  increasing  criticality  of  these 
systems  to  Orders  of  Battle  (OOBs)  demands  an  increased 
appreciation  of  the  effects  of  the  ionosphere  on  electro- 
magnetic signals-  Additionally,  there  is  growing  awareness 
that  attention  should  be  f ocussed ,  not  only  on  the  well- 
covered  and  well -understood  middle  latitudes,  but  also  on 
potential  conflict  areas  in  the  high  latitudes,  such  as  the 
so-called  Greenland- Iceland-United  Kingdom  (G-I-UK)  Gap  and 
the  Northern  Pacific-  Unfortunately,  the  high  latitude 
ionosphere  (>65°)  is  probably  the  most  complex,  and  in  the 
polar  cap  in  particular,  the  least  understood  and  least 
studied  of  the  ionospheric  domains.  Recent  discoveries  of 
anomalies  in  the  polar  cap,  as  well  as  the  growing 
sophistication  and  availability  of  data  collection 
techniques,  has  contributed  to  a  surge  of  interest  in  the 
high  latitude  ionosphere. 

Of  particular  interest  are  the  anomalies  known  as  polar 
cap  arcs.  Previous  studies  of  the  arcs  have  identified 
their  structure  and  dynamics  (Weber  and  Buchau ,  1981); 
energy  characteristics  of  the  electron  population  (Hardy, 
1984);   Interplanetary   Magnetic    Field    (IMF)   dependency 


(Gussenhoven, 1982) ;  and  various  characteristics  of  the 
magnetosphere- ionosphere  coupling  processes  related  to  the 
arcs. 

The  Air  Force  Geophysics  Laboratory  (AFGL) ,  among 
others,  is  conducting  ongoing  investigations  in  an  effort  to 
understand  and  characterize  the  high  latitude  ionosphere. 
Data  analyzed  in  this  thesis  was  collected  from  the  AFGL 
stations  in  Thule  Air  Base  and  Qanaaq ,  Greenland,  during  the 
months  of  January  and  February,  1985.  The  data  features 
recordings  of  high  elevation  angle  satellite  signals  and 
photometer  images  taken  simultaneously.  The  intent  of  the 
analysis  in  this  thesis  is  to  discover  if  the  appearance  of 
polar  cap  arcs  is  associated  with  any  changes  in  the  total 
electron  content  (TEC)  or  any  amplitude  scintillation 
problems  in  the  satellite  signal. 


1 1 .   BACKGROUND 


The  relationship  between  the  ionosphere  and  electro- 
magnetic radiation  is  a  particularly  intimate  one:  solar 
radiation  in  the  extreme  ultraviolet  (EUV)  and  X-ray  range 
creates  and  maintains  the  ionosphere,  while  earth-based 
radio  frequency  (RF)  radiation  determined  its  existence, 
studies  it  and  uses  it  as  a  vital  transmission  medium. 

Although  the  existence  of  an  electrically  conducting 
layer  of  the  atmosphere  had  been  postulated  earlier, 
confirmatory  evidence  became  possible  only  with  the 
development  of  radio,  beginning  around  the  turn  of  the 
century.  Since  then,  RF  has  provided  an  important  means  of 
studying  the  ionosphere,  beginning  with  the  relatively 
simple,  but  effective  and  still  widely  employed,  technique 
of  using  the  high  frequency  (HF)  band  (~2-30  MHZ)  for 
vertical  sounding  to  determine  the  boundaries  of  the 
ionosphere.  Increasingly  sophisticated  techniques  are 
associated  with  progressively  higher  bands.  For  example, 
incoherent  scatter  radar,  used  to  determine  electron 
temperature  and  density,  employs  the  very  high  frequency 
(VHF)  and  ultra  high  frequency  (UHF)  bands  (~30-300  MHz  and 
"300  MHz-3  GHz);  and  the  analysis  of  perturbations  of 
satellite  signals  in  the  upper  end  of  the  UHF  range  (L-band 
microwave)  is  used  to  determine  total  electron  content. 
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A.  MORPHOLOGY  OF  THE  IONOSPHERE 

The  ionosphere  is  a  blanket  of  stratified,  partially 
ionized  plasma  surrounding  the  earth  at  an  altitude  of  "50- 
1000*  kilometers.  Ionospheric  behavior  is  governed  by 
processes  which  result  in  the  production/destruction  and 
transport  of  ionization;  different  mechanisms  dominate 
according  to  both  altitude  and  latitude.  "Layering"  into 
the  D,  E,  Fl  and  F2  regions  occurs  because  (a)  the  ionizing 
energy  of  solar  radiation--a  production  mechanism- -is 
deposited  according  to  the  absorption  characteristics  of  the 
atmospheric  constituents;  (b)  recombination- -a  destruction 
mechanism- -depends  on  atmospheric  density;  and  (c)  the 
composition  of  the  atmosphere  changes  with  altitude,  and 
prevailing  types  of  production  and  loss,  as  well  as 
transport,  mechanisms  are  tied  to  these  changes  in 
composition.  Prevailing  mechanisms  also  change  with  respect 
to  the  geomagnetic  latitude,  reflecting  the  influence  of  the 
earth's  magnetic  field. 

B.  PRODUCTION  AND  LOSS  IN  THE  D,  E,  AND  Fl  REGIONS 

The  principal  production  process  in  the  lower  and  middle 
ionosphere  at  low  and  middle  latitudes  (<65°)  is  photo- 
ionization.  Photo- ionization  occurs  when  energy  from  solar 
EUV  and  X-ray  photons  is  absorbed  by  atmospheric 
constituents,  resulting  in  the  creation  of  ion-electron 
pairs.    Since   both   the   density   and   composition   of  the 
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atmosphere  change   with  altitude,  different  constituents  are 
ionized  by  different  wavelengths  and  energies  of  radiation. 

The  D  region  ("50-90  km)  is  characterized  by  positive 
ions  of  nitric  oxide  (NO*  )  and  molecular  oxygen  (02*  ),  with 
Lyman  a  and  hard  X-rays  being  the  main  solar  source  of 
ionizing  radiation.  Molecular  nitrogen  (N2)  is  also  ionized 
in  significant  amounts,  but  is  rapidly  converted  to  02*  . 

The  E  region  ("90-140  km)  is  characterized  by  02*  and 
N0+,  with  some  contribution  from  N2*  and  atomic  oxygen  (0*  ). 
The  solar  radiation  bands  responsible  are  chiefly  Lyman  0, 
the  Lyman  continuum,  Carbon  III,  and  soft  X-rays. 

In  the  F  region  ( > 140  km),  0*  is  the  predominant  ion, 
but  also  present  are  02*  ,  N2*  ,  He*  ,  and  N*  .  The  Lyman 
continuum  and  radiation  at  <800  3.  are  the  principal  sources 
of  ionizing  energy. 

Photo- ionization  is  still  the  more  prolific  ion 
producing  process  at  high  latitudes,  but  in  this  instance, 
corpuscular  ionization  makes  a  significant  contribution  to 
the  ion  population,  particularly  at  night,  when  it  is  the 
only  source  of  ionization.  Corpuscular  ionization  results 
from  the  interaction  (i.e.,  collisions)  of  precipitating 
charged  particles  and  the  atmospheric  constituents. 
Particle  precipitation  occurs  toward  the  geomagnetic  poles, 
where  the  geomagnetic  field  lines  are  closest  to  vertical. 

The  "destruction"  of  ionizat ion--loss  mechanisms--dif fer 
particularly  with  altitude.   In  the  D  region,  the  lowest  and 
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most  dense,  a  number  of  processes  seem  to  be  operative.  Of 
these,  the  most  important  are  believed  to  be  dissociative 
recombination  of  electrons  with  hydrated  water  cluster  ions 
and  photodetachment .  Loss  processes  in  the  E  and  Fl  region 
are  better  known  and  are  clearly  dominated  by  dissociatve 
recombination  of  02*  and  NO* .  Table  1  summarizes  the 
characteristics  of  the  layers  (Chamberlain,  1978,  p. 165). 

C.   THE  F2  LAYER 

The  F2  is  the  highest  of  the  distinguishable  ionospheric 
regions.  The  dominant  production  process  is  still  photo- 
ionization.  Because  the  region  is  optically  thin, 
attenuation  is  not  a  factor  as  in  the  lower  regions,  and  the 
ionization  rate  of  the  constituents  is  proportional  to  their 
number  densities.  In  the  F2  layer,  the  major  constituent 
ionized  is   0,  with  a  relatively  small  contribution  from  02. 

Above  a  certain  altitude,  known  as  the  F2  peak,  neutral 
density  is  low  enough- -and  with  electron  density  decreasing 
--that  plasma  diffusion  processes,  rather  than  photochemical 
loss  mechanisms  characterize  the  behavior  of  the  region.  In 
the  lower  F2  layer,  below  the  F2  peak,  both  electron  and 
neutral  densities  are  high  enough  to  provide  a  basis  for 
photochemical  loss  mechanisms. 

The  controlling  loss  process  is  the  very  rapid  ion-atom 
or  ion-molecule  charge  exchange  reaction  (which  converts  0* 
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to  NO*  and  02*),  followed  by  dissociative  recombination. 
Because  the  rate  of  the  reaction  is  constrained  by  the 
neutral  density,  which  decreases  with  altitude,  the 
occurrence  of  the  loss  process  decreases  exponentially  with 
height,  and  electron  density  continues  to  increase  beyond 
the  altitude  of  peak  production- 
After  a  certain  height,  the  plasma  diffusion  rate 
becomes  more  rapid  than  the  photochemical  loss  rate,  and 
dominates  the  upper  F2  layer,  causing  the  electron  density 
to  decrease  with  height. 

The  F2  peak  is  the  height  at  which  the  electron  density 
is  highest.  It  is  a  result  of  the  net  effect  of  the  loss 
process  and  the  plasma  diffusion  process.  The  peak  occurs 
at  the  level  where  the  time  constants  for  the  plasma 
diffusion  and  photochemical  loss  rates  are  approximately 
equal-   (Gnanalingam,  1987) 

D.   THE  HIGH  LATITUDE  IONOSPHERE 

The  high  latitude  plasma  population  is  primarily 
governed  by  solar  EUV  ionization  and  precipitating 
electrons,  as  well  as  horizontal  and  magnetic  field  aligned 
transport  processes. 

The  solar  EUV  component  varies  smoothly.  However, 
because  the  magnetosphere- ionosphere  (M-I)  coupling 
processes  are  both  dynamic  and  interactive,  electron 
precipitation  throughout  the  high  latitude  ionosphere  varies 
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greatly  in  energy  distribution  and  with  latitude  and  local 
time.  Field  aligned  currents,  electric  fields  and  charged 
particles  represent  the  principal  coupling  mechanisms 
between  the  magnetosphere  and  the  ionosphere.  Electric 
fields  are  responsible  for  the  large  scale  plasma  convection 
patterns  that  are  characteristic  of  the  high  latitude 
regions.  Charged  particles,  precipitating  along  magnetic 
field  lines  from  the  magnetosphere,  are  believed  to  be  the 
source  of  the  upward  field  aligned  currents,  which  are  also 
associated  with  gradients  in  the  electric  fields.  The 
particles  may  also  be  accelerated;  they  typically  have  a 
higher  range  of  energies  than  that  of  their  solar 
wind/magnetosheath  source  region.  This  results  in  not  only 
horizontal  and  vertical  plasma  structures  but  also  a  source 
of  free  energy  for  the  plasma  instability  processes 
characteristic  of  the  region. 

The  major  domains  of  the  high  latitude  ionosphere- -the 
auroral  zone  and  the  polar  cap--are  generally  identified  by 
both  geomagnetic  latitude  and  characteristic  phenomena. 
They  differ  in  magnetic  field  topology,  convective  electric 
field  patterns  and  particle   precipitation  characteristics. 

The  auroral  zone  encompasses  the  region  from  ~60°-75° 
Corrected  Geomagnetic  Latitude  (CGL).  The  magnetic  field 
lines  are,  for  the  most  part,  closed,  and  the  source  of  the 
precipitating  electrons  is  the  magnetosphere ' s  plasma  sheet: 
the  central   plasma  sheet  (CPS)  for  the  equatorward  portion, 
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with  particle  energies  up  to  a  few  kilo  electron  volts 
(keV);  and  the  boundary  plasma  sheet  (BPS)  for  the  poleward 
portion,  with  particle  energies  up  to  a  few  hundred  eV .  The 
large  scale  electric  fields  are  set  up  to  support  sunward 
convection,  except  in  the  Harang  discontinuity  region  near 
midnight  and  the  dayside  cusp  near  noon. 

The  polar  cap  is  the  poleward  region  bounded  by  the 
instantaneous  auroral  oval.  Magnetic  field  lines  are  open, 
or  close  at  such  great  distances  down  the  magnetotail  that 
they  may  be  treated  as  open.  Polar  electron  precipitation 
is  believed  to  originate  directly  from  the  magnetosheath , 
and  is  customarily  classed  by  energy  level.  Polar  rain  has 
typical  energies  of  "100  eV  and  is  uniformly  distributed 
over  large  areas.  Polar  showers  have  mean  energies  of  wl 
keV ,  and  are  enhanced  flux  areas  embedded  in  larger  regions 
of  polar  rain.  Polar  squalls  are  localized  zones  of  intense 
precipitation,  with  particle  energies  of  up  to  several  keV . 
The  polar  cap  is  generally  an  area  of  ant i -sunward 
convection.   (Burke,  1982,  pp. 685-686;  Heelis,  1982,  p.567) 

E.   POLAR  CAP  IRREGULARITIES 

The  high  latitude  F  layer  is  riddled  with  irregularities 
which  have  a  broad  range  of  scale  sizes  (from  10"  2  -  103  m) 
and  a  global  distribution.  Interest  in  the  high  latitude 
ionospheric  domain  has  lately  begun  to  focus  on  the 
heretofore  neglected  polar  cap  as  well  as  the  more  obviously 
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dynamic  auroral  zone.  Recently  recognized  irregularities  in 
the  polar  cap  illustrate  both  particle  precipitation  and  the 
convective  transport  processes  that  dominate  the  region. 

F- layer  patches  are  regions  of  enhanced  ionization 
produced  in  the  dayside  ionosphere,  which  are  then 
transported  into  and  across  the  polar  cap  and  ultimately 
into  the  nightside  auroral  oval  (Figure  lb).  The  patches 
are  characterized  by  airglow  emission  at  6300  &;  high 
electron  density  (5-10  times  background  level);  large  size 
(~1000km  in  horizontal  extent);  and  speeds  of  up  to  1000 
m/sec  in  an  anti-sunward  direction.  The  patches  also  have 
an  IMF  dependency,  occurring  when  the  IMF  has  a  southerly 
component  (i.e.,  Bz<0)  and  when  the  geomagnetic  field  is 
moderately  disturbed.  As  these  patches  move  through  the 
ionosphere,  they  cause  significant  increases  in  the  ambient 
TEC  and  measurable  amplitude  scintillation  in  both  UHF  and 
L-band  satellite  signals  through  which  they  pass.  (Weber  et 
al.,  1986,  p. 12121) 

A  second  recently  recognized  phenomenon  is  the  polar  cap 
arc  (Figure  2).  The  arcs  tend  to  appear  during  magnetically 
quiet  periods,  and  when  the  IMF  has  a  northerly  component 
(i.e.,  Bz >0) .  They  are  nominally  sun-aligned,  and  drift  in 
a  direction  perpendicular  to  their  alignment,  that  is,  dawn- 
to-dusk  (Figure  la).  They  have  a  curtain-like  structure: 
height  and  length  are  large  relative  to  width.   They  stretch 
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from  an  altitude  of  ~250kra  to  "1000* km.  As  in  the  patches, 
the  arcs  have  an  optical  signature  at  6300  & ,  and  are 
associated  with  upward  field  aligned  currents.  (Akasofu, 
1985,  pp. 15-18;  Weber  and  Buchau,  1981,  p. 125;  Carlson  et 
al.,  1984,  p. 895  and  p. 898) 

F.   DESCRIPTION  OF  THE  EXPERIMENT 

Data  used  in  this  thesis  was  collected  at  the  AFGL 
stations  at  Thule  and  Qanaaq,  Greenland.  The  instruments 
used  were  an  All  Sky  Imaging  Photometer  (ASIP),  at  Qanaaq, 
and  a  specially  equipped  Global  Positioning  System  (GPS) 
satellite  receiver  at  Thule.  Thule  is  located  at  76°32'N, 
68°42'W,  and  Qanaaq  at  77°28'N,  69°13'W  (geographic 
coordinates).  Their  great  circle  separation  is  104.5  km, 
and  the  true  azimuth  angle  from  Qanaaq  to  Thule  is  179°39'. 
Although  the  ASIP  and  the  GPS  receiver  are  not  collocated, 
their  spatial  separation  is  not  considered  significant 
relative  to  the  large  area  covered  by  the  ASIP  (Figure  3) . 

Polar  cap  arcs  are  believed  to  be  caused  by 
precipitating  soft  electrons,  and  would  therefore  represent 
regions  of  increased  electron  density  and  TEC.  Ionospheric 
sounders  are  the  most  common  means  for  remotely  determining 
electron  density.  However,  the  arcs  cause  spread  F  on 
sounder  data,  making  it  extremely  difficult  to  extract 
correct  critical  frequency  and  thus  determine  electron 
density.   An  entirely  different  technique  is  used  in  the 
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Figure  3.  Geometry  of  GPS  receiver  and  ASIP  locations 


collection   of   GPS   data,   and   using   the  GPS  data  allows 
unambiguous  determination  of  TEC. 

The  spread  F  masking  in  the  ionogram  indicates  the 
presence  of  structured  irregularities  in  the  arcs. 
Irregularities  also  cause  amplitude  scintillation,  and  the 
occurrence  of  changes  in  scintillation  indices  caused  by  the 
arcs  are  further  indication  of  the  presence  of 
irregularities.  The  GPS  data  includes  amplitude 
scintillation  measurements,  taken  simultaneously  with  the 
TEC  measurements,  and  allows  further  definition  of  the  arcs 
and  the  irregularities  they  contain. 

At  the  time  the  images  and  satellite  receiver  recording 
were  made  (January  and  February),  the  Thule/Qanaaq  area  was 
undergoing  24  hour  periods  of  darkness.  Therefore,  solar 
EUV   was    minimal.      The   ASIP   is   programmed   to   image 
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continuously  under  moonless  conditions  after  the  sun  has 
reached  a  solar  depression  angle  of  10°,  and  data  from  the 
satellite  receiver  is  recorded  continuously.  Of  the  5  GPS 
satellites  currently  operating,  at  least  1  was  in  view  up  to 
4  hours  at  an  elevation  angle  of  greater  than  20° ,  for  a 
total  of  ~18  hours. 

G.   THE  ALL  SKY  IMAGING  PHOTOMETER 

The  ASIP  is  a  monochromatic  imaging  system  with  high 
sensitivity  (threshold  =  ~20  Rayleighs)  developed  for  low 
light  level  research.  An  upward- looking,  wide  angle  lens 
images  a  large  section  of  the  sky  (155°  field  of  view)  onto 
an  interference  filter  near  the  focal  plane.  The  optics  are 
telecentric  in  design,  ensuring  the  same  size  light  cone  at 
the  focal  plane  for  each  point  in  the  field  of  view,  and 
allowing  the  use  of  narrowband  interference  filters  ("25  2). 
Although  up  to  four  filters  can  be  mounted  on  a  filter  wheel 
to  allow  sequential  measurements,  only  the  6300  &  01  and  the 
4278  3>  M2*  were  used  at  Qanaaq.  The  field  is  re- imaged  onto 
the  faceplate  of  a  two  stage  image  intensifier.  The  image 
is  then  coupled  to  a  35  mm  Automax  camera,  using  Kodak  4X 
film  at  "10  seconds  exposure. 

Filters  at  6300  &  and  4278  A  are  used  because  those 
wavelengths  tag  the  photochemical  reactions  of  specific 
atoms  and  molecules  which  predominate  at  different  heights 
in  the   ionosphere,  and  can  therefore  identify  the  altitudes 
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at  which  the  phenomena  occur.  For  example,  arcs  have  been 
observed  at  both  E  and  F  region  heights,  and  their  presence 
can  be  optically  verified  because  the  filters  monitor  the 
following   reactions: 

soft  particles   >  e"  f  *  *  *  +  0 >0*  ~«~-»  hi?  6300  a 

(which  occurs  at  F  region  heights  where  the  dominant 
atmospheric  constituent  is  0) ,  and 

harder  particles  -->  e"'iSt  +  N2  >  N2*  ~— »  h^  4278  a 

(which  occurs  at  E  region  heights  where  N2  is  the  dominant 
atmospheric  constituent). 

Control  sequencing  is  provided  by  a  programmable  camera 
control  system  which  controls  the  mechanical  operation  and 
camera  exposure  setting.  Date,  time  and  filter  setting  are 
automatically  recorded  on  each  image.  The  ASIP  also  has  a 
light  sensor  which  will  disable  the  camera  under  moonlit  or 
sunlit  conditions. 

The  ASIP  is  most  often  used  to  map  the  structure  and 
dynamics  of  auroral  and  airglow  structures  rather  than  to 
measure  quantitatively  absolute  intensities.  In  this  case, 
the  ASIP  was  used  to  verify  the  presence  of  an  arc  and  to 
indicate  its  position  in  relation  to  the  ray  path  of  the 
satellite.  (Weber  et  al . ,  1977, pp. 7-9) 
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Figure   4   lists   system   specifications   (Weber  et  al., 
1977,  p. 11). 


Field  of  View: 

155° 

Pass  Band: 

25  A  at  f  1.4;  5  A  at  f  8 

Resolution: 

1/2     zenith,    2     horizon 

Spectral  Response: 

S  -  20,    exceeding  100  ^A /lumen 

Picture  Storage: 

No  detectable  degradation  for  up  to  3  sec 

Tube  Gain: 

Photon  noise  granularity  visible  above  tube  noise 

Threshold  Sensitivity: 

20  R  at  2  sec  exposure 

Dynamic  Range: 

20  R  to  10  kR  covered  by  3  preset  HV  settings 

Flatness  of  Field: 

30  percent  loss  at  edge  of  field 

Reoetition  Rate: 

Tvoicallv  20  sec  for  complete  filter  cycle 

Figure  4.   Summary  of  ASIP  specifications 

H.    THE  GPS  RECEIVER  AND  MEASUREMENT  TECHNIQUES 

The  NAVSTAR  Global  Positioning  System  (GPS)  is  a 
satellite  based  navigation  system  which  provides  accurate  3- 
dimehsional  position  and  timing  information  to  users. 

The  constellation,  eventually  planned  to  contain  18 
satellites,  currently  has  5,  at  63°  inclination  in  12  hour, 
sidereal -synchronous  orbits.  The  satellites  broadcast 
continuously  on  1.57542  GHz  (LI)  and  1.22760  GHz  (L2). 
Their  signals  are  coherently  derived,  biphase  modulated  at 
10  MHz  with  a  pseudo- random  noise  (PRN)  code  and  contain 
ephemeris,  satellite  clock  drift,  signal  code  state  and 
ionospheric  propagation  delay  data.  The  user's  receiver 
measures  the  pseudo  range  (measured  transit  time  of  the 
signal  multiplied   by  the   speed  of   light)  and  pseudo  range 
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rate  ,  using  the  encoded  signals  from  4  satellites.  These 
measurements  can  then  be  converted  to  a  user  3-dimensional 
position  and  velocity.  (Demaro,  1981,  pp. 36-37;  Cretcher, 
1975,  pp. 406 -407) 

The  advantages  of  the  GPS  satellites  for  high  latitude 
research  lie  in  (a)  their  inclination,  which  allows  high 
look  angles  for  significant  periods  of  time  in  the  true 
polar  regions,  well  above  the  cusp  and  the  auroral  oval,  and 
(b)  the  system's  2  frequency  design  feature  which  allows  the 
user  to  correct  for  signal  delays  due  to  ionospheric 
effects.  (Demaro,  1981,  pp. 36-37;  Cretcher,  1975,  pp. 406- 
407)  . 

One  of  the  critical  ionospheric  effects  on  radio  wave 
propagation  is  the  group  delay  of  the  modulation  envelope  of 
the  wave-  Other  ionospheric  effects,  measurable  by  the  2 
frequency  design  of  the  GPS  system,  are  the  differential 
carrier  phase  advance  and  amplitude  and  phase  scintillation. 
The  ionospheric  parameter  responsible  for  these  effects  is 
TEC,  the  total  number  of  free  electrons  which  the 
transmitted  signal  encounters  in  the  ionosphere.  The  group 
delay  and  carrier  phase  advance  are  proportional,  to  the 
first  order,  to  the  TEC  in  the  transionospheric  path  of  the 
wave.   (Klobuchar,  1983,  pp. 5  and  16) 
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1 .   Group  delay 

Group  delay  is  the  additional  time  delay  of  a 
signal  transmitted  from  above  the  ionosphere,  and  is 
generally  given  by: 

At  =  (40.3/cf2  )  x  TEC 
where  c   is  the   speed  of   light,  f   is  the  frequency  of  the 
wave,  and  TEC  is  in  el/m2 . 

For  a  system  using  a  pair  of  fairly  widely  spaced 
frequencies,  group  delay  as  a  function  of  frequency  is: 
At,  =  (k/cf,  2)  x  TEC  ,  dt2  =  (k/cf22)  x  TEC 
where  t,  is  the  group  delay  at  frequency  f,  and  t2  is  the 
group  delay  at  frequency  f 2 .  Iff,  is  the  primary  frequency 
and  f2  is  the  lower  frequency  chosen  for  correction  of 
ionosphe-rically  induced  errors,  then: 

S(At)  =  (k/c)  x  TEC  (l/f22  -  1/f,  2) 
=  At,  (f,  2  -  f22  )/f22 
or 

At,  =  (f22/f,  2  -  f22)  x  S  (At) 

where  S(At)  is  the   difference  in   the  time   delays  measured 
simultaneously  on  the  2  frequencies,  in  nanoseconds. 

With  GPS,  the  user's  receiver  converts  the 
received  spectrum  to  equivalent  carriers  with  10.23  MHz 
modulation.  The  10  MHz  modulation  is  transmitted  with  a 
known  phase  difference  on  the  2  carriers,  and  the  received 
modulation  phase   difference  is   a  direct  measurement  of  the 
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group  delay.    The   group  delay  is  related  to  the  equivalent 
TEC  by: 

TEC  =  (2.852  x  10' 6  )  x  5  ( At ) 
To  obtain  the  absolute  TEC,  one  value  of  TEC  is 
obtained  to  fix  the  relative  phase  advance  scale  to  an 
absolute  one  throughout  each  satellite  pass.  Usually, 
several  independent  values  of  group  delay  are  fitted  on  an 
r.m.s.  basis  to  the  continuous  relative  phase  advance  scale. 
(Klobuchar,  1983,  pp. 7-9;  Klobuchar  et  al.,  1985,  p. 2) 
2.   Differential  carrier  phase  advance 

In  addition  to  group  (or  modulation)  delay,  the 
phase  of  the  carrier  is  changed  by  the  ionosphere,  with  the 
carrier  phase  being  advanced  with  respect  to  its  phase  in 
the  absence  of  an  ionosphere.  Carrier  phase  advance  is 
related  to  TEC  by: 

5  C+) )  =  (1.34  X  lCT"/f)  x  TEC 
where  f  is  the  carrier  frequency  in  Hertz. 

The  GPS  satellites  transmit  2  coherently  derived 
carrier  frequencies  for  the  purpose  of  measuring 
differential  carrier  phase  advance.  For  the  GPS  frequencies 
LI  (1.575  GHz)  and  L2  (1.228  GHz),  the  phase  advance 
difference,  referenced  to  the  lower  frequency,  is: 

A(Scp)  =  (1.34  x  107/L2)  x  (m2-l/m2)  x  TEC 
where   m   is   L1/L2   for   the   GPS  frequencies.   For  the  GPS 
f requenc  ies , 

^(Scp)  =  (4.31  x  10"  '  7  )  x  TEC 
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or 

TEC  =  (2.32  x  10' 6  el/m2 )  Z*CS<p) 
(per  complete  2tt  cycle  of  differential  carrier  phase  advance 
between  LI  and  L2,  measured  at  L2)  - 

Because  the  2  coherent  GPS  carriers  are  relatively 
close,  the  measured  ionospheric  phase  advance  at  L2  minus 
that  at  LI  is  only  39%  of  the  total  phase  advance  at  L2. 
Still,  the  technique  has  a  sensitivity  of  better  than  0.1 
radians  r.m.s.  in  a  16  Hz  bandwidth  using  an  omnidirectional 
antenna,  and  has  the  added  advantage  of  removing 
contributions  to  the  phase  advance  due  to  geometric  Doppler 
changes.  (Klobuchar,  1983,  p. 10;  Klobuchar  et  aL,  1985,  p. 2) 

The  differential  carrier  phase  advance  measurement 
accurately  measures  relative  changes  in  TEC-  Background 
TEC,  however,  generally  produces  much  more  than  one  complete 
cycle  of  carrier  phase  advance.  Group  delay  is  used, 
therefore,  to  measure  absolute  TEC. 
3 .   S« 

The  third  parameter  measured  by  the  GPS  receiver  is 
amplitude  scintillation,  which  is  presented  as  changes  in 
the  S«  index.  The  S,  scintillation  index  is  the  normalized 
variance  of  signal  intensity  and  is  a  widely  used  measure  of 
signal  fluctuations  caused  by  ionospheric  irregularit ies. 
(Yeh  and  Liu,  1982,  p. 324) 
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I.   REDUCTION  OF  GPS  DATA 

All  the  GPS  data  were  recorded  on  5  channel  magnetic 
tape  and  processed  at  AFGL  using  a  computer  program 
developed  specifically  for  that  purpose.  The  program 
unpacks  and  reformats  the  data  for  compatability  with  the 
AFGL  CYBER  computer,  determines  calibration  levels  for  each 
channel,  converts  voltages  to  TEC  units,  and  produces  3 
output  files  for  additional  analysis.  The  first  of  these 
contains  6  second  averages  of  group  delay  (absolute  TEC). 
The  second  file  contains  S«  values  measured  at  L2  at  1 
minute  sliding  averages.  The  third  file  contains  deramped 
carrier  phase  advance  values  (relative  TEC).  (Doherty,  1987) 
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III.   METHODOLOGY  OF  THE  EXPERIMENT 

Data  simultaneously  recorded  by  the  ASIP  and  the  GPS 
receiver  yielded  four  types  of  information:  (a)  6300  £> 
optical  images;  (b)  absolute  TEC;  (c)  relative  TEC;  and  (d) 
S* .  ASIP  images  were  used  to  determine  both  the  presence  of 
the  arcs  and  the  periods  when  they  crossed  the  raypath  of 
the  GPS  receiver.  Data  from  the  GPS  receiver  were  then 
examined  to  establish  the  degree  of  correlation,  if  any, 
between  the  presence  of  an  arc  and  changes  in  the  TEC  levels 
and  the  S«  index. 

A.   ANALYSIS  OF  DATA 

Comparison  of  data  showed  excellent  correlation  between 
the  onset  of  an  arc's  passage  through  the  GPS  raypath  (i.e., 
the  arc's  leading  edge)  and  a  relative  increase  in  TEC.  A 
total  of  7  arcs  from  January  and  February  1985  were 
examined;  4  are  represented  here.  For  the  instances 
examined,  the  average  relative  increase  in  TEC  was  ~3.6 
units  (1  unit  =  101  6  el/m2),  with  the  lowest  increase  at 
"1.7  units  and  the  greatest  at  ~5.0  units.  In  all  but  one 
of  the  instances,  there  was  also  a  small  but  measurable 
change  in  the  S«  index:  an  average  increase  in  S«  from  ~0.06 
to  ^0.1  (  S*  =  ~0.04).  Further  comparison  showed,  however, 
that  the  arcs  followed  different  patterns   of  behavior,  even 
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with  respect  to  the  clear  association  between  onset  of 
leading  edge  and  TEC  enhancement.  Specific  cases  will  be 
presented  to  illustrate  some  of  the  differing 
characteristics  of  the  arcs. 

1.   Case  1 :  26  January  1985 

The  primary  arc  of  this  event  began  encounter  with 
the  GPS  signal  at  ~0804  UT  on  26  January  1985,  while  still 
relatively  unformed.  The  arc  cohered  in  the  raypath,  and 
did  not  clear  it  until  0920  UT.  The  images  (Figures  5  and 
6)  show  the  slow  formation  and  stagnation  of  the  arc  and 
this  is  reflected  in  the  clear  cut  but  gradual  increase  in 
TEC.  The  peak  relative  increase  in  TEC  was  ^4.4  units, 
which  occurred  at  0822  UT,  early  in  the  arc's  transit.  The 
arc  hovered  in  or  near  the  raypath  for  over  an  hour,  but  the 
TEC  levels,  while  remaining  elevated,  never  reached  the 
relative  maximum  achieved  at  0822  UT.  A  plot  of  the 
relative  TEC  changes  is  at  Figure  7. 

In  this  instance,  the  end  of  encounter  is  also 
unmistakable.  The  arc  had  completely  cleared  the  raypath  by 
0920  UT,  and  the  relative  TEC  levels  dropped  abruptly. 
Although  not  shown  here,  the  relative  TEC  levels  continued 
undisturbed  for  some  hours  after  the  arc's  passage. 

Figure  8  is  a  plot  of  the  relative  TEC  for  the 
same  25  period  on  the  following  day  (0729  -  0947  UT,  27 
January  1985) ,  during  which  there  were  no  arcs  recorded  or 
observed. 
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The  GPS  satellites  repeat  their  orbits  with  only  a  four 
minute  daily  change  in  their  rise  and  set  times,  so  that 
Figures  7  and  8  represent  data  taken  from  the  same  satellite 
at  very  nearly  the  same  times  on  subsequent  days. 
Comparison  of  the  TEC  curves  clearly  reveals  significant 
differences  in  the  TEC  behavior  on  the  two  days.  During  the 
passage  of  the  arc  on  26  January,  the  TEC  variations  are 
many,  large  and  abrupt,  reflecting  the  disturbed  conditions 
generated  by  the  arc.  The  27  January  plot  shows  only  a  very 
gradual  and  very  small  decrease  from  0804  -  0920  UT ,  with  no 
enhancements  to  mar  the  smoothness  of  the  curve.  Since  this 
tallies  with  the  fact  that  there  were  no  arcs  observed,  this 
plot  is  illustrative  of  conditions  unperturbed  by  the 
presence  or  passage  of  an  arc. 

The  obvious  increase  in  TEC  was  not,  in  this 
instance,  coupled  to  any  measurable  change  in  S,  (Figures  9 
and  10). 

2.   Case  2:  16  February  1985 

This  arc  also  had  distinct  leading  and  trailing 
edges,  but  its  transit  was  much  faster  than  that  of  the  arc 
just  examined-  It  entered  the  raypath  already  well  formed, 
at  1037  UT,  16  February,  had  traversed  it  completely  and 
dissipated  (out  of  the  raypath)  by  1053  UT.  It  has  a  less 
substantial  appearance  in  the  ASIP  images  (Figure  11),  and 
maximum  relative  change  in  TEC  was  "2.2  units  (Figure  12). 
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Although  the  change  in  TEC  was  relatively  small, 
the  arc's  transit  did  cause  an  appreciable  change  in  the  S«  , 
an  increase  from  ~0.06  to  ~0.12  (Figure  13). 

An  unusual  feature  of  this  arc  was  that  its 
lifetime  in  the  field  of  view  was  as  a  single  arc.  In 
general,  the  arcs  seem  to  arrive  or  form  parallel  to  each 
other  in  groups,  which  gives  a  streaky  or  striated 
appearance  to  the  ASIP  image.  Alternatively,  a  single  arc 
will  divide  itself  to  form  two  or  more  separate  arcs. 
3.   Case  3:  25  January  and  16  February  1985 

There  were  many  instances  which  exhibited  the 
distinct  leading  edge  but  which  have  highly  ambiguous  end 
conditions.  These  arcs  either  have  no  distinguishable 
trailing  edge,  or  are  followed  so  closely  by  other  arcs  or 
turbulence  that  they  can  be  classed  as  one  event. 

The  first  example  of  this  began  at  2330  UT,  25 
January,  as  indicated  the  extremely  rapid  rise  in  TEC,  with 
a  peak  increase  in  level  of  "4.7  units-  The  ASIP  image 
shows  a  discernable  arc  moving  through  the  raypath  until 
"2340  UT,  when  the  primary  arc  had  passed.  Almost 
immediately,  however,  the  images  reveal  the  raypath  embedded 
in  a  secondary  arc  which  began  to  form  out  of  the  main  arc 
at  "2338  UT  (Figure  14).  This  secondary,  which  does  not 
eventually  form  a  discernable  arc,  remained  around  the 
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raypath  for  a  number  of  minutes  before  moving  on,  and  caused 
a  series  of  relative  enhancements  of  TEC,  with  the  maximum 
relative  increase  at  ^5.0  units  (Figure  15).  The  primary 
arc  caused  a  noticable  increase  in  S«  ,  from  ~0.6  to  ^0.16, 
while  the  disturbed  area  following  it  generated  an  S« 
increase  of  "0.04  (Figure  16). 

A  similar  but  not  as  severe  instance  began  at  0911 
UT  on  16  February  (Figure  17).  The  22  minute  transit  of  the 
primary  arc  caused  relative  TEC  enhancements  with  a  maximum 
value  of  "3.8  units.  Again,  however,  the  primary  arc  is 
closely  followed  by  its  own  partially  formed,  split  off 
secondary  and  an  area  of  generalized  turbulence  that  never 
successfully  coalesces.  Although  the  relative  TEC 
enhancement  of  the  follow-on  area  is  modest,  possibly 
because  true  minimum  was  not  reached  between  the  arcs,  the 
TEC  is  increased  to  a  level  which  is  equal  to  that  caused  by 
the  primary  arc  (Figure  18).  With  the  gradual  dissipation 
of  the  follow-on  area,  the  TEC  declined  to  pre-event  levels. 

The  S«  index  for  this  event  was  also  disturbed 
(Figure  19).  The  disturbances  were  small  but  still  obvious: 
a  0.03  increase  for  the  primary  arc  and  an  increase  of  ~0.02 
for  its  follow-on  area. 
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Figure  17.  ASIP  Images:  0851  -  0959  UT ,  16  February  1985 
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IV.   DISCUSSION 


Arcs  and  patches  are  the  only  structured  anomalies  to 
have  been  discovered  in  the  polar  cap.  Since  the  arcs  occur 
when  B*>0,  and  since  Bz >0  conditions  exist  approximately  50% 
of  the  time,  the  arcs  can  be  said  to  account  for  roughly 
half  the  irregularity  incidences  in  the  polar  cap. 
Therefore,  if  the  arcs  are  associated  with,  or  indicative 
of,  ionospheric  irregularities,  their  occurrence  is  a 
significant  factor  in  calculating  and  predicting  ionospheric 
effects  on  C3 I  systems. 

Two  results  seem  conclusive  from  the  analysis  of  data 
done  for  this  thesis.  First,  the  arcs  are  clearly 
associated  with  large  and  rapid  variation  in  ambient  TEC. 
Absolute  TEC  values  were  not  calibrated  for  the  period 
during  which  these  arc  observations  were  made.  However, 
they  were  calibrated  from  GPS  code  data  for  a  similar 
observation  period  a  year  earlier  (January  -  February  1984), 
when  solar  flux  was  nearly  the  same.  During  magnetically 
quiet  periods,  the  absolute  equivalent  vertical  TEC  values 
were  less  than  10  x  10' s  el/m2  ,  as  shown  in  Figure  20. 
Figure  20  (Klobuchar  et  al . ,  1986,  Figure  3)  compares 
theoretical  calculations  of  polar  cap  TEC  (indicated  on  the 
graph  by  the  solid  dots),  which  take  into  account  solar  EUV , 
average   polar   rain   precipitation,   and   the   earth's   co- 
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rotation;  and  actual  recorded  values  of  TEC  for  the  quiet 
day  of  1  February  1984.  If  the  absolute  TEC  values  given  in 
Figure  20  are  assumed  to  be  approximately  the  same  as  for 
the  period  January  -  February  1985,  then  the  observed  TEC 
increases  associated  with  the  arcs  can  be  as  much  as  100%. 
(Klobuchar  et  al . ,  1986;  Klobuchar,  1987) 
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Figure  20.   Theoretical  (solid  dots)  and  measured   absolute 
equivalent  vertical  TEC 

Second,  the  arcs  also  produce  small  but  measurable 
amplitude  scintillation  effects  even  at  the  comparatively 
high  GPS  L2  frequency.  The  frequency  and  intensity  of 
ionospheric  irregularities  increase  with  sunspot  number,  so 
amplitude  scintillation  effects  vary  with  the  level  of  solar 
activity.  These  arcs  were  observed  during  a  period  close  to 
solar  minimum  (which  occurred  in  1986);  S«  variations  are 
expected  to  be  much  greater  in  the  years  approaching  solar 
maximum. 

Table  2  is  a  summary  of  the  TEC  and  S«  changes  in  the 
arcs  discussed. 
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V.   SUMMARY 


First-order  effects  of  ionospheric  irregularities 
decrease  as  frequency  increases.  Generally,  then,  HF 
systems  are  more  likely  than  L-band  systems  to  experience 
directly  ionospherically  induced  outages.  Even  in  the  VHF 
range,  however,  where  the  majority  of  U.S.  satellite 
communications  frequencies  are  located,  irregularity  caused 
effects,  specifically  amplitude  scintillation,  are  bad  and 
frequent.  Deep  signal  fades  of  >20dB  are  common  and  well 
documented  (Basu  et  al . ,  1986,  p. 20  and  Figure  10). 
Additionally,  as  frequencies  are  increased  to  accomodate 
satellite  based  systems,  particularly  in  communications, 
flexibility  in  the  system  is  lost  almost  completely. 

As  frequency  goes  up,  the  accuracy  threshold  for 
minimum  acceptable  performance  is  also  raised.  L-band 
signals  will  not  be  greatly  degraded  by  amplitude 
scintillation,  for  example,  but  the  selectivity  of  the 
receivers  must  be  so  narrow  to  deliver  accurate  processed 
information  that  even  small  perturbations  may  cause  a 
degraded  product.  A  system  such  as  GPS,  which  advertises 
extreme  accuracy  and  incorporates  self -correct  ion  for 
ionospheric  time  delay  as  a  basis  of  the  system,  is 
particularly  vulnerable  to  this  kind  of  error. 
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Also  susceptible  are  the  space-based  radar  systems, 
particularly  synthetic  aperture  radar  (SAR).  SAR  requires  a 
high  degree  of  phase  coherence  for  imaging  integrity.  The 
phase  coherence  of  a  SAR  signal  can  be  easily  distorted  by 
phase  scintillation  effects  introduced  by  irregularities 
with  kilometer  range  scale  sizes,  such  as  those  associated 
with  polar  cap  arcs.  The  coherency  break-up  of  a  SAR  signal 
cannot  be  easily  compensated  for,  if  at  all,  in  the  signal 
processing,  since  irregularities  are  neither  spatially  nor 
temporally  predictable.  ( Szuszczewicz  et  al . ,  1983,  pp. 766- 
768) 

For  the  bulk  of  the  C3 I  systems  in  use  or  under 
development,  then,  the  behavior  of  the  F  layer  is  an 
important- -sometimes  essential- -factor  in  determining  how 
well  those  systems  perform.  With  the  increasing  strategic 
and  tactical  importance  of  the  very  high  latitude  areas,  the 
significance  of  the  behavior  of  polar  cap  F  layer  is  also 
increased. 

F  layer  polar  cap  arcs,  which  generate  enhancements  in 
TEC  and  measurable  amplitude  scintillation  effects, 
therefore  have  the  potential  for  degrading,  sometimes 
seriously,  the  performance  of  C3 I  systems  operating  in  or 
across  the  polar  latitudes. 
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